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Fabien Jean/’^ Tao Zhou,^’^ Nils Blanc,Roberto Felici,^ Johann Coraux,^’^ and Gilles Renaud^’^ 

^CNRS, Inst NEEL, F-38042 Grenoble, France 
^ Univ. Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France 
^GEA, INAG-SP2M, Grenoble, F-38054, France 
‘^European Synchrotron Radiation Facility, Boite Postale 220, F-38043 Grenoble Gedex 9, France 

(Dated: June 25, 2015) 

The structure of a graphene monolayer on Ir(lll) has been investigated in situ in the growth 
chamber by surface x-ray diffraction including the specular rod, which allows disentangling the 
effect of the sample roughness from that of the nanorippling of graphene and iridium along the 
moir-like pattern between graphene and Ir(lll). Accordingly we are able to provide precise esti¬ 
mates of the undulation associated with this nanorippling, which is small in this weakly interacting 
graphene/metal system and thus proved difficult to assess in the past. The nanoripplings of graphene 
and iridium are found in phase, i.e. the in-plane position of their height maxima coincide, but the 
amplitude of the height modulation is much larger for graphene (0.379 ±0.044 A) than, e.g., for the 
topmost Ir layer (0.017 ± 0.002 A). The average graphene-Ir distance is found to be 3.38 ± 0.04 A. 


Graphene, a monoatomic layer of carbon atoms 
arranged in a honeycomb lattice, has been investi¬ 
gated thoroughly in the past ten years because of 
its exceptional properties, which hold promises for 
numerous applicationsii Transition metal surfaces 
form a broad family of substrates for the growth 
of large area, high quality graphenei^ New proper¬ 
ties can be induced in graphene through the interac¬ 
tion with the substrate, e.g. electronic bandgaps^ 
spin-polarization^ and superconductivity.— In most 
graphene-on-metal systems, the interaction is mod¬ 
ulated at the nanoscale, due to lattice mismatch be¬ 
tween graphene and the metal, which results in two- 
dimensional patterns with periodicity of the order 
of nanometers, often referred to as ’’moirs”, follow¬ 
ing an analogy with the beating of optical waves 
through two mismatched periodic lattices {e.g. tis¬ 
sue veils). Knowledge on the topographic properties 
of these moirs, i.e. the average graphene-metal dis¬ 
tance, and the perpendicular-to-the-surface ampli¬ 
tude of the graphene and metal undulations across 
the moir, is desirable in view of characterizing the 
interaction and rationalizing the other properties. 

The topography is however hard to grasp at 
such small scales. Most efforts that have relied on 
scanning tunnelling microscopy have faced the is¬ 
sue of the entanglement of the structural and lo¬ 
cal density of state which is inherent to the tun¬ 
nelling effect. A striking illustration has been 
the debate on the amplitude^ii and signSiS of the 
moir-related undulation in graphene/Ru(0001) and 
graphene/Ir(lll), respectively. Atomic force mi¬ 
croscopy has proven to be a valuable alternative, 
provided it is performed with care, especially with 
respect to the possible chemical interaction between 



FIG. 1. Sketch of the reciprocal space, the hexagonal 
grid shows the partition of its {II,K) plane according 
to the lO-on-9 commensurability. H, K and L are in 
reciprocal lattice unit of the moir (superlattice) surface 
unit cell. In gray are shown the measured CTRs from 
the iridium, with circles to highlight the positions of the 
different Bragg reflections. The graphene rods are shown 
in black. The specular CTR (H—K—O) is shown in red. 
Each is labeled with its {H,K) position in the lO-on-9 
moir surface supercell. 


tip and sample— Siii Scattering techniques, such as 
low-energy electron diffraction (FEED), surface X- 
ray diffraction (SXRD), and X-ray standing waves 
(XSW), are free of such probe-induced perturba¬ 
tions of the systems. To the expense of com¬ 
plex calculations in the framework of the dynam¬ 
ical theory of diffraction, LEED was used to as¬ 
sess the topography of graphene/Ru(0001)^ and 
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FIG. 2. Experimental structure factors Fh,k{L) of irid¬ 
ium CTRs and graphene rods from SXRD measurements 
of the first sample in black with the error bars. The solid 
red lines represent the best fit with the Fourier model. In 
blue with the rods is the contribution of a flat graphene 
layer alone, to highlight the effect of the roughness and 
undulations on the rods. The specular rod (0,0) from the 
second sample is reported in the bottom right in black 
with error bars. The solid red line represents the final fit, 
the green one the contribution of the iridium alone and 
in blue the contribution of a flat graphene layer alone. 


graphene/Ir(lll)fii SXRD was used to analyze the 
topography of graphene/Ru(0001)fi^ as was done by 
XSW for graphene/Ir(lll).— Confirming and refin¬ 
ing the results obtained with these approaches is 
of prime importance in order to set reliable points 
of reference for first principle calculations, which 
are cumbersome in essence in such systems due to 
the importance of non-local {e.g. van der Waals) 
interactions 

Here, we address the model graphene/Ir(III) sys¬ 
tem, typical of a weak graphene-metal interaction. 
Its moir topography only slightly deviates from the 
flat case and is thus difficult to characterize. With 
the help of two techniques, SXRD and extended x- 
ray reflectivity (EXRR), the latter not having been 
employed to characterize monolayer graphene on a 
substrate before, we deduce an average 3.38±0.04 A 


distance between graphene and Ir(lll), and deter¬ 
mine, with an uncertainty as low as with scanning 
probe microscopyA^ a 0.379 ± 0.044 A amplitude of 
the graphene undulation. Besides, we are able to es¬ 
timate the undulation of the Ir layers, which is usu¬ 
ally not accessible to other techniques, 0.017 ±0.002 
A for the topmost one. 

The synchrotron x-ray diffraction measurements 
were performed in ultra-high vacuum chambers cou¬ 
pled with z axis diffractometers at the BM32 and 
ID03 beamlines of the European Synchrotron Ra¬ 
diation Eacility. Details on the chambers and the 
beam are given in Ref. [isl . The non-specular crys¬ 
tal truncation rods (CTRs) were measured on BM32 
and the specular rod, OOL, was measured by EXRR 
on ID03. The x-ray beam energy was set at II keV. 
The reciprocal space scans of the scattered inten¬ 
sity presented below are all normalized to the in¬ 
tensity measured with a monitor placed before the 
sample. Eor the SXRD measurements, the intensity 
along the Ir(III) crystal truncation rods (CTRs) 
and along the graphene rods was measured with 
a Maxipix two-dimensional detector in stationary 
mode for the upper range of the out-of-plane scat¬ 
tering vector component {i.e. large values of the 
out-of-plane reciprocal space coordinate L), and by 
performing sample rocking scans for low I^valuesA^ 
The amplitude of the structure factors Fh,k{L) - the 
square root of the measured intensity - for the dif¬ 
ferent CTRs and graphene rods, corresponding each 
to different values of the in-plane reciprocal space 
parameters H and K, were extracted and processed 
with the PyRod program described in Ref. [l^. Py- 
Rod was also used to simulate the structure factors 
using the model described below, and to refine the 
structural parameters of this model with the help of 
a least squares fit of the simulation to the data. The 
total uncertainty on the experimental structure fac¬ 
tors is dominated by the systematic error estimated 
to be 6.1%, according to Ref. [2 the statistical error 
being everywhere smaller than I %. 

The Ir single crystals were cleaned according to a 
procedure described in Ref. [l^ allowing for consid¬ 
erably reducing the concentration of residual carbon 
in bulk Ir(lll). Graphene was grown in two steps, 
first by 1473 K annealing of a room-temperature- 
adsorbed monolayer of ethylene, second by exposure 
to 10“® mbar of ethylene the surface held at 1273 
K. This growth procedure allows for selecting a well- 
defined crystallographic orientation of purely single¬ 
layer graphene with respect to IrGll).— Compared 
to the samples studied in Refs. 0 and m the sur¬ 
face coverage is larger (100%) in our case. Our 
growth procedure is similar to that used to prepare 
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the 100%-coverage graphene studied in Ref. 0, yet 
the temperatures which we chose for each step are 
different, actually identical to those used for prepar¬ 
ing one of the samples addressed in Ref. E We note 
that both the graphene coverage and growth temper¬ 
ature have been argued to influence the structure of 
graphene^*^ and thus its propertiesi^ Two samples 
were prepared, one in each of the UHV chambers in¬ 
stalled at the BM32 and ID03 beamlines where the 
SXRD and EXRR experiments were performed re¬ 
spectively. The hexagonal lattice unit cell of the 
iridium surface has a lattice parameter of 2.7147 A 
at room temperature. The graphene unit cell has a 
measured lattice parameter of 2.4530 A. The ratio 
between the two lattice parameters, 0.903, is close 
to 0.9. Therefore, in the following we assume that 
the system is commensurate, with a (1010) graphene 
cell coinciding with a (99) iridium one. In the fol¬ 
lowing, the in-plane unit cell of reciprocal space is 
the moir onei^ This corresponds to 77 or if indexes 
multiples of 9 and 10 for Ir CTRs and graphene rods, 
respectively (Fig. [T|). 

Figure [2] shows the Ir CTRs and graphene rods. 
As expected for a (essentially) two-dimensional layer 
such as graphene, the graphene rods are basically 
featurelessj^ Qualitatively, because the undulations 
of the graphene and top substrate layers are ex¬ 
pected to be small, the main features are i) the 
pronounced interference effect on the specular rod 
Fo,o{L) related to the average distance dzQr between 
Ir and graphene, expected to be larger than the bulk 
distance of 2.2 A; ii) the decrease of the otherwise 
featureless CTRs in between Bragg peaks, related 
to the substrate roughness ; and Hi) the decrease of 
the graphene rods with increasing T, dominated by 
the undulation of the graphene layer, as shown with 
the simulated graphene rods for a flat graphene layer 
alone in Fig. [5] This decorrelation between rough¬ 
ness and undulation allows these parameters to be 
determined with high accuracy. 

In order to achieve a quantitative characteriza¬ 
tion of the topography of the system, we introduce 
a simple model<^ A limited set of parameters (see 
Fig. m, including the average interplanar distances, 
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the actual roughness, and the amplitude of undu¬ 
lation of each layer, seems to be a reasonable op¬ 
tion for a simple modeling of the system. In order 
to approach this description, we introduce a lattice 
model based on a Fourier series, such as the one pro¬ 
posed for graphene/Ru(0001)i^ In this model, the 
displacement in the direction i {i = {x,y,z\) of an 
atom with x, y and z coordinates, with respect to 
the corresponding position in a flat layer, is given 
Gr ; 

Ir, ; 

Ir2 ^ 

Irs ■ 

0---Bulk Ir 

FIG. 3. Sketch of the parameters studied. In black is 
the graphene and in blue, red and green are the surface 
layers of iridium. The amplitudes of their corrugation 
are shown by arrows in the middle. The start of the 
bulk iridium is sketched with the dotted black line. The 
gray dashed lines represents the expected bulk positions 
for the different atomic layers without corrugation. The 
z-axis on the left is a reference to the linear dependency 
of the iridium corrugation amplitude. 
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by 

dr* = A* ^xsin[27r(sa;-|-t^)]-|-ii* ^xcos[27r(sx-|-i^)] 

S,t 

( 1 ) 

where the sum runs over the different orders of 
the series. Due to the crystal symmetry of graphene 
and Ir(lll), the displacements must respect a p3ml 
symmetry, i.e. they must fulfill 


Rj-Hdr[Rj(r)]} = Rj{dr[Rj-yr)]} (2) 

with j G [0,5]. Rq is the identity matrix, Ri and 
R 2 correspond to the ± 120 rotations and the last 
three to the mirror planes. 


, R.3 





( 3 ) 


Fourier coefficients in Eq. (1) are independent. 


These symmetry constraints impose that not all 
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Their relationships are given in Table HI 

In the following we further simplify the model 
by limiting the Fourier development to first order, 
which is legitimate due to the fact that no significant 
diffraction data is measurable beyond first order (a 
diffraction experiment is actually a measurement of 
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the Fourier transform of the electronic density, thus, 
to a good approximation, of the shape of graphene). 
Besides, we assume that the undulations of all Ir 
layers are in phase, as found in density functionnal 
theory (DFT) calculations. In this framework, the 
X, y and z displacements simply write: 


dr'^ = X (2 X sin(27ra:) + sin(27r2/) + sin(27r(x — y))) 

+ X (2 X cos(27rx) — cos(27r?/) — cos(27r(a; — y))) 


dr^ = X (sin(27ra:) + 2 x sin(27r7/) — sin(27r(a; — y))) 

+ X {cos{2'kx) — 2 X cos(27ry) + cos(27r(a: — y))) 


( 5 ) 


dr^ = X {2 X sin(27ra;) — 2 x sin(27r?/) — 2 x sin(27r(a; — y))) 

+ X (2 X cos(27ra;) + 2 x cos(27ry) + 2 x cos(27r(x — y))) 


( 6 ) 


Thus, only two variables per atomic plane, A^ and 
, are needed to describe the in-plane displace¬ 
ments. The model is applied to graphene/Ir(III), 
with three iridium layers and one graphene layer. 
Each of these layers is characterized by four Fourier 
coefficients {A^, B^, A^ and B*), plus another pa¬ 
rameter corresponding to an average z displacement 
of the layer from its equilibrium position in the bulk. 
This distance between metal planes parallel to the 
surface is known to vary, in some cases by as much 
as few percents and in a non monotonous manner 
across the few topmost layers of metal surfaces^^. In 
order to reduce the number of free parameters how¬ 
ever, we assume a linear dependence of the distance 
between Ir(lll) planes, thus of of A^ and B^), as a 
function of depth. This assumption complies with 
the results of the DFT simulations (cf. Table [TTl 
where the topographic parameters of the model are 
listed). 

The Fourier model was used to fit the SXRD data. 
The expected in-plane displacements (Fig.|4]), below 
O.OI A according to first principle calculations?^ have 
no noticeable influence on the Ir CTRs and graphene 
rods, and are discarded in the simulations!^ The 
best fit lead to a value of 3.5 and the results 
are shown in Table |lTl We find a 98 ± 2% graphene 
coverage. The graphene is found to have a mean 
distance of dzcr = 3.39 ± 0.28 A with its substrate 
and a corrugation of Azor = 0.379 ± 0.044 A. The 
graphene distance with its substrate is close to the 
interlayer spacing in graphite, 3.36 A. As explained 
above, the benefit of the SXRD analysis of both 


graphene and Ir contributions is to provide an ac¬ 
curate value of the amplitude of the graphene un¬ 
dulation perpendicular to the surface, as compared 
to other techniques. The interlayer Ir spacings are 
found to be 2.203 ± 0.012 A, 2.212 ± 0.007 A and 
2.223 ± 0.002 A from top to bottom. The topmost 
layer of iridium has an undulation of 0.017±0.002 A, 
the second layer has an undulation of O.OlliO.OOl 
A, and the last one is 0.006 ± 0.001 A. Finally, the 
roughness of the iridium substrate is found to be 
0.42 ± 0.20 A, following a simple /3-modeli^ This 
small value may be linked with the small coherence 
length of the X-ray beam (corresponding to about 
10 flat Ir terraces separated by atomic step edges) 
on the BM32 beamline. 

atop fee hep atop 
0.379 .. - .. 

0.017AI MMMMMMMMMMS 
0.011 AI MMMMMMMMMMS 

0.006 AI •0mJKm090909Cm090K» 


FIG. 4. Cut of the lO-on-9 commensurability to repre¬ 
sent the corrugations and displacements of the atomic 
layers. Carbon atoms of the graphene are black circles, 
the iridium atoms are in blue, red and green to show the 
ABC stacking of the different layers. The three coinci¬ 
dence regions of graphene with the substrate as well as 
the corrugations and interlayer spacing are denoted. 


3.38 A 
2.203 A 
2.205 A 
2.225 A 


The best fit between simulations and SXRD data 
is achieved for an iridium undulation in phase with 
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the graphene one, with a smaller amplitude though. 
This finding is at variance with that obtained in ear¬ 
lier scanning probe microscopy measurements per¬ 
formed in specific imaging conditions?^ and supports 
the picture progressively assembled through other 
reports, based on scanning probe microscopies 
XSW?^ and first principle calculations 

The main limitation of this SXRD analysis is 
the rather large uncertainty on the dzcr distance. 
This motivated complementary measurements of the 
specular rod on the second sample, using the ID03 
setup as the extended reflectivity was not accessible 
in the BM32 setup. The EXRR result is shown in 
Fig. [5] together with the best fit and simulated and 
graphene specular rods. The best fit of the specular 
rod, yielding a value of 1.064, was done with a 
simplified model, in which the undulations of both 
the iridium or graphene were fixed at the values ob¬ 
tained from the SXRD analysis. It yields a value 
of dzQr = 3.38 A, very close to that determined on 
the other sampler by off-specular SXRD, but with a 
much better accuracy, ± 0.04 A. The graphene layer 
of this second sample is found incomplete, with a 90 
± 2% graphene coverage. In addition, the spacings 
between the topmost Ir planes, found to be 2.203 
± 0.010 A, 2.205 ± 0.008 A and 2.225 ± 0.004 A 
from top to bottom (2.217 A in the bulk). We find 
interlayer distances between Ir(lll) planes that are 
larger than those found in the absence of graphene 
This finding is consistent with the p-doping found 
for graphene^ which implies electron transfers from 
graphene to Ir(lll). The corresponding higher elec¬ 
tronic density in Ir(lll) is expected to counterbal¬ 
ance the surface relaxation in bare Ir(lll). The sub¬ 
strate roughness in this case is found to be 1.1 ± 0.1 
A, larger than that obtained from SXRD. This is 
however expected since the coherence length of the 
beam is two orders of magnitude larger here, e.g. 
around 1000 atomic steps of the substrate scatter 
the beam coherently. 

This is the first study of a sample with a complete 
graphene coverage, thus the deviations from previ¬ 
ous studies can be explained due to strains in the full 
layer that can relaxe in graphene island. This could 
also be explained by the difference in the growth 
process (temperature, methods...) and Busse et al^ 
showed that the undulation varies depending on the 
graphene coverage. Moreover, the undulation could 
also be affected by the growth methods (full/partial 
growth, chemical vapor deposition, temperature pro¬ 
grammed growth...) and growth temperature as it 
has been reported that these parameters affect the 
graphene lattice parameter and its commensurabil- 
ity with the substrateA^d^ The iridium undulations 


are also found larger than those deduced from a 
FEED studyiii This might be due to some limita¬ 
tion of LEED to analyse layers below the graphene 
one, because of the small electron mean free path. 

The graphene-metal distance which we obtain is 
in good agreement with values deduced by XSW, 
LEED, and AFM (see Table HIl) . The undulation of 
the graphene which we obtain is also in agreement 
with that found by LEED and AFM. It is however 
smaller than that deduced from XSW. The differ¬ 
ence might originate from two effects. First, we 
recently found that the in-plane lattice parameter 
of graphene varies as a function of the preparation 
method, which is different in Refs. 1^ andfllL several 
TPG cycles at 1420 K for different coverages and 
one TPG at 1500 K respectively, and in the present 
work, TPG at 1473 K followed by CVD at 1273 K for 
complete coverage. Given that the strain is closely 
related to the graphene buckling (undulation) we 
indeed expect different undulations in each of these 
reports. Second, the strain (and thus buckling) of 
graphene was argued to depend on the fraction of 
edge atoms in graphene, i.e. on gr^hene coverage^ 
Our results, unlike those in Refs. 0 and [H address 
close-to-full layer graphene. 

The Fourier model was also tested to fit the dis¬ 
placements obtained by the DFT calculations de¬ 
scribed in Ref. The model was in very good accor¬ 
dance with the DFT calculations results, in particu¬ 
lar the iridium top layer and graphene, thus confirm¬ 
ing that the first order Fourier component is enough 
to describe the system, as shown in Fig [5] More¬ 
over, it also confirmed that and of the three 
iridium layers have an almost linear dependence as 
a function of depth. From the DFT simulation, the 
corrugations of the iridium surface layers from top 
to bottom are 0.015 A, 0.012 A and 0.04 A while 
the graphene one is 0.35 A, which are close to the 
experimental results. 

In fact, this analysis has a limite too, as our start¬ 
ing hypothesis on the structure of the supercell, a 
(1010) graphene cell coinciding with a (99) irid¬ 
ium, may have an impact on the results. It was 
reported previously that this system cannot be con¬ 
sider fully commensurate, as it is really a compo¬ 
sition of commensurate domains with incommensu¬ 
rate boundaries'^ and that the thermal history of the 
sample effects itA^d^ Here, the 9.03 ratio indicates 
that there should be a combinaison of (I0I0)/(99), 
(2I2I)/(I9I9) and incommensurate domains. How¬ 
ever, despite the complexity of the sample, the start¬ 
ing hypothesis of the problem allows to extract a 
good approximation of the actual structure. 

To conclude we have employed SXRD to deter- 
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mine with high resolution, on the basis of a simple 
structural model, the structure of a weakly scatter¬ 
ing atomically thin membrane, graphene, in weak 
interaction with a metallic substrate made of strong 
scatterers, Ir atoms. We determine the undula¬ 
tion of graphene across the moir-like superstructure 
formed between graphene and Ir(lll), 0.379 ±0.044 
without the ambiguity inherent to other ensemble¬ 
averaging techniques. Our determination of the av¬ 
erage graphene-Ir(lll) distance is consistent with 
previous reports based on local-probe and ensemble¬ 
averaging analysis. Finally we unveil the faint cor¬ 
rugations predicted by DFT calculation in the sub¬ 
strate, which are as low as 0.017 ± 0.002 for the 


topmost Ir layer, and are characteristic of a weak 
C-Ir bonding having a slight covalent character in 
some of the sites of the moir. The use of SXRD for 
other two-dimensional membranes, such as transi¬ 
tion metal dichalcogenides, boron nitride, or mono- 
layer silica, should allow for constructing a compre¬ 
hensive picture of the nanomechanics of atomically 
thin membranes under the influence of substrates. 

We thank Olivier Geaymond, Thomas Dufrane 
and the staff members of the ID03 and BM32 
beamlines, Nicolae Atodiresei for the DFT calcu¬ 
lations data and the French Agence Nationale de 
la Recherche for funding (Contract No. ANR-2010- 
BLAN-1019-NMGEM). 
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TABLE I. Relationships between the Fourier coefficients Alt and Bit {i = {a, y,z}) 



SXRD (I'** sample) 

EXRR (2"‘‘ sample) 

DFT 

Ref. Ig 

Ref. [ill 

dZGr 

3.39 ±0.28 

3.38 ± 0.04 

3.43 

3.38 ± 0.04 

3.39 ±0.03 

AzOr 

0.379 ± 0.044 


0.46 

0.6 ±0.1 

0.47 ±0.05 





1.0 ±0.2 


dZiri 

2.203 ±0.012 

2.203 ±0.010 

2.190 


2.222 

dZlr2 

2.212 ±0.007 

2.205 ± 0.008 

2.175 


2.224 

dZirs 

2.223 ± 0.002 

2.225 ± 0.004 

2.184 


2.222 


0.017 ±0.002 


0.015 


0.006 


0.011 ±0.001 


0.012 


0.006 

^^ZXt‘ 

0.006 ± 0.001 


0.004 


0 

p 

0.42 ± 0.20 

1.05 ±0.08 




Ocr 

98 ± 2% 

89.7 ± 1% 

100% 

39% 

Partial 





63% 



TABLE 11. Topographic parameters for the two samples, the DFT calculations data from Ref. and results from 
Ref. (XSW) andii (LEED + AFM). dzcr is the mean distance between the graphene and its substrate; Azcr is 
the graphene undulation amplitude; dzin, dzir 2 dzir^ are the interlayer distances of the iridium surface layers 
and Aziri 1 Azin ^^d Azir-^^ are their undulation amplitudes; p is the roughness of the sample surface; Oor is the 
graphene coverage in percent. All the parameters are in ngstrms (A) except the coverage. 



o 

Q. 

N 



c 
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Q. 

N 


Fourier series model experimental fit 




FIG. 5. Sketches of the graphene out of plane variations from (a) the Fourier model used on the experimental data 
and (b) the DFT calculations results and the out of plane variations of the topmost iridium layer from (c) the Fourier 
model used on the experimental data and (d) the DFT calculations results, with the DFT calculations results shown 
in left half-discs and the Fourier series fit in right half-discs The out of plane corrugation is shown with a color 
gradient, with the scales in A. 


































